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ABSTRACT 

We present a three-pointing study of the molecular gas in the starburst nucleus of M82 based on 190- 
307 GHz spectra obtained with Z-Spec at the Caltech Submillimeter Observatory. We present intensity 
measurements, detections and upper limits, for 20 transitions, including several new detections of CS, 
HNC, C2H, H2CO, and CH3CCH lines. We combine our measurements with previously-published 
measurements at other frequencies for HCN, HNC, CS, C"^^S, and HCO+ in a multi-species likelihood 
analysis constraining gas mass, density and temperature, and the species' relative abundances. We 
find some 1.7-2.7x 10* of gas with riHa between l-6xl0'*cm~^ and T >50 K. While the mass and 
temperature are comparable to values inferred from mid- J CO transitions, the thermal pressure is a 
factor of 10-20 greater. The molecular interstellar medium is largely fragmented and is subject to 
ultraviolet irradiation from the star clusters. It is also likely subject to cosmic rays and mechanical 
energy input from the supernovae, and is warmer on average than the molecular gas in the massive 
star formation regions in the Milky Way. The typical conditions in the dense gas in M82's central 
kpc appear unfavorable for further star formation; if any appreciable stellar populations are currently 
forming, they are likely biased against low mass stars, producing a top-heavy initial mass function. 
Subject headings: galaxies: abundances — galaxies: individual (M82) — galaxies: ISM — galaxies: 
starburst — instrumentation: spectrographs — techniques: spectroscopic 



1. INTRODUCTION 

Studies of molecular gas - the actively star-forming 
part of the interstellar medium (ISM) - in other galaxies 
have been carried out most extensively in the low-lying 
rotational transitions of the CO molecule. Due to its high 
abundance relative to molecular hydrogen and its favor- 
able energy level spacing as compared to typical molecu- 
lar cloud temperatures, CO produces the brightest lines 
in the millimeter-wavelength portion of the spectrum. 

However, while these lines trace the bulk of the 
molecular gas, the small dipole moment of CO re- 
sults in modest critical densities for these lines, ~ 
10"^ — 10^ cm~^, and thus they do not strongly dis- 
criminate between high-density star-forming cores and 
more diffuse gas. High-dipole-moment molecules, such 
as HCN and CS, have critical densities and radiative 
rates that are 100 — 1000 times larger than CO. De- 
spite their much smaller abundances, these species can 
be used to probe the dense (lO'^ — lO'' cm~^) cloud 
cores believed to be associated with star formation (SF). 
Early mea surement of HCN and CS in the Galaxy 
(IMorris et'al. 1 9 74j: | Turner et aIlll973D aiid in external 
galaxies (jHenkel et al .lll988l:ISolomon et al.ll99l showed 
that their intensities are well-correlated with the total 
far-infrared (FIR) flux. More recently, HCN J = 1 ^ 
luminosity has been shown to be directly proportional to 
FIR luminosity (a proxy for star formation rate [SFR]) 
in a sample of ^100 normal spiral and starburst galax- 
ies ([Gao fc Solomon. 2004ai .b). as well as individual star- 
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formation sites in the Galaxy (jWu et al.l 12005'). 

Insofar as HCN J = 1 — >■ measures dense gas mass, 
the LpiR / HCN correlation across 7-8 orders of magni- 
tude in luminosity implies a scale-independent relation- 
ship between dense gas mass and SFR. An accurate as- 
sessment of the physical conditions in the HCN-emitting 
gas as well as the mass scaling (AfH2/-^HCN) is thus of 
universal interest for theoretical SF studies. Is the HCN- 
emitting gas simply a bi-product of star formation, per- 
haps formed in shocks or outflows and not participating 
in the formation of new stars, or is some of it the very 
material from which new stars form? Such an assess- 
ment is best made with multiple transitions of HCN and 
by including transitions from other high-dipole-moment 
molecules where available. 

To assess the average properties of the dense gas on 
the scale of a nuclear starburst, we have observed the nu- 
cleus of the M82 in a suite of millimeter-wave transitions 
of high-dipole moment species. The brightest infrared 
(IR) galaxy in the sk y due to its proximity (3.9 Mpc, 
iSakai fc Madorel[l999f) . M82 radiates an inf rared lumi- 
nosity (L = 5.9xlO^°L0. ISanders et"ani2003[ ). exceeding 
that of the Galaxy, from a r egion that is only about 450 
pc in radius (e.g., Leeuw & Robson 2009). Because of 
this concentration of star-forming activity, M82 has been 
dubbed the prototypical starburst galaxy, which makes 
it a particularly interesting laboratory for the study of 
SF. It has been suggested for over 30 years that the stel- 
lar initial mass function (IMF) in M82 (and presumably 
other starburst nuclei) is biased ag ainst low -mass stars 
relative to the the loca l IMF ( Rieke et al.l [T98a ,19931 
iForster Schreiber et al.|[2 003). bu t this ha.s been debated 
(e.g.. ISatvapal et al.l 11997; .ColbeTt et al.lll999D . If the 
IMF is indeed low-mass deficient, a plausible line of in- 
quiry is the initial conditions of SF - the temperature 
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and density of the dense molecular cloud cores. 

We have obtained full 190-307 GHz spectra with the 
Z-Spec 1-mm grating spectrometer, which accesses the 
J = 3 ^ 2 transitions of HCN, HCO+, and HNC and 
J = 4->3, J = 5-^4, and J = 6 -> 5 transi- 
tions of CS, with a uniform calibration. Our study 
benefits from the prior observat ions of J = 1 — »■ 
trans i tions of HCN, HNC HCO+ (|Nguven-Q-Rieu et al.l 
119891: iHuettemeister et al.lll995D. as well as J = 4 -» 3 
transitions of HCN, HCO+ (|Seaauist fc FraveH l200l . 
These data are combined with our observations to gen- 
erate the first comprehensive multi-species excitation 
and radiative transfer model for the dense gas in this 
source. Of course, since M82 is also wel l- studied in 
muhiple CO transitions f Wild et al.l '1992^; 'Ma o et al 



pOO; Ward et al. 2003; Wc ii3 et al.l 12005: SeaquisleLal 
[2006; Pan uzzo et al.l 12010) as well as in the mid- and 



far-IR atomic gas t racers (jForster Schreiber et al.l [20031 : 
IColbert et an[l999l) . we have the opportunity to put the 
dense gas into context with the other ISM components, 
as well as the general properties of this prototypical nu- 
clear starburst. 

2. OBSERVATIONS 

Z-Spec is the first grating spectrometer for the mil- 
limeter band; it covers the full 190-307 GHz range 
instantaneously, dispersing this band to an array of 

16 bolometer s. Mo re information can be found 

in iGlenn et al.l ((200l . [Bradford etall ([20091 ) and the 
Society of Photo- Optical Instrument a tion Engineers 
(SPIE) articles: ([Navlor et al.l [200l [Bradford et all 
I2004t lEarle et al.ll2006( ) The instrument operates at the 
Nasmyth focus of the Caltech Submillimeter Observa- 
tory (CSO) atop Mauna Kea. While the instrument is 
1// stable down to ^ 100 mHz, we use a chop and nod 
mode to avoid the atmospheric fluctuations, which be- 
come important relative to the fundamental noise sources 
at ^0.3-1 Hz, depending on the atmospheric conditions. 
For the M82 observations the chop frequency was 1.6 Hz, 
the throw was 90" in azimuth, and the nod interval was 
20 seconds. Three pointings along M82's major axis 
were observed on 2009 January 5 as summarized in Fig- 
ure [T] and Table [TJ The M82 spectra are calibrated using 
Mars, with an interpolation scheme using bolomet er op- 
erating voltages as a measure of response ( Earl j [20081 : 
[Bradford et all [2001 . The data were reduced with stan- 
dard demodulation and differencing appropriate for the 
chop and nod observing mode and Uranus is used as a 
spectral flat-fielder. We expect the channel-to-channel 
calibration uncertainties to be less that 10% except at 
the lowest frequencies which are extremely sensitive to 
the wing of the 186 GHz atmospheric water line. 

3. RESULTS 

Spectra for the NE, GEN, and SW pointings are shown 
in Figure [2] The general features agree with previous 
measurements of the 1.2-mm continuum and CO J = 
2^-1, such as those in Figure [T] 

3.1. Continuum Fluxes and Fits 

The continuum fiux in the Z-Spec band is due to a com- 
bination of the thermal dust emission which dominates 
at higher frequencies and the free-free emission from [H 



Table 1 

M82 Observed Positions 



Pointing 


a offset 
(arcsec) 


5 offset 
(arcsec) 


Int. Time 
(min) 


Sensitivity 
(Jy sV2) 


NE 


-1-12.2 


-1-3.3 


60.2 


1.3 


CEN 


-1-2.7 


-0.5 


68.6 


1.4 


SW 


-6.1 


-3.8 


60.1 


1.2 



Note. — The R. A. and Dec. offsets are relative to 
OJ2000.0 = 9''55'"51.9=, <5j2000.o = 69°40'47.14". The 
integration time is the total demodulated time; the on- 
source time is half the listed values. The quoted sensitiv- 
ity is the median value of channel errors multiplied by the 
square root of the integration time and does not represent 
the ultimate sensitivity of the instrument; Z-Spec's sensi- 
tivity to spectral lines and to fainter sources is better by 
at least a factor of two. The optical depth during these 
observations was T225GHz = 0.08 — 0.1. 



II] region electrons which dominates at lower frequencies. 
The dust emission is parameterized as 



Friiy) ^ nBy{T) <^ 1 - cxp 



(1) 



where B JT) = (2hv^ /f )/icx-p\hiy/kT]-l) is the Planck 
function. iHughes et al.l y,994) fit a collection of observa- 
tions from 3.3 mm to 40 /zm and derive fl = 1.34 x 10~* 
sr, T = 48.1 K, /3 = 1.3 and Aq = 7.9 microns for the 
thermal dust. The free- free emission is given by a simple 
power law, 



(2) 



Fs{i^)=Fo — 



where i^o is the free- free fiux at i/q. Analysis of 
a map of M82 at 92 GHz found that the emission 
at that frequency is domina ted by free-free emission. 
iCarlstrom &: Kronbergl ([19911 ) calculate the free-free flux 
density Fff(92 GHz)= 0.5 ± 0.1 Jy. In the Z-Spec band, 
the free-free emission accounts for roughly one quarter of 
the total continuum flux, but that is signiflcant enough 
that a dust-only fit to our spectra does not match the 
continuum level. Instead, we model the continuum flux 
we observe to be a beam-scaled fraction of the total con- 
tinuum emission of the galaxy. 



F{v) = A 



240 GHz 



[Ft{v) + Fs{y)] 



(3) 



This model assumes that the free-free and thermal con- 
tinuum emission have the same spatial distribution, 
which is reasonable given our coarse spatial resolution. It 
also accounts for the spectral dependence caused by the 
spatial distribution. If M82 completely fllled our beam 
at all frequencies, we would expect B = while B = 2 
would be the prediction for a point source. We find that 
i? = 1 is the best fit value for all three pointings (see 
Table [2]) , which makes sense given that M82 is observed 
edge on and is thus roughly point-like in one dimension 
and beam-filling in the other, relative to our 30" beam. 
The continuum fraction seen in the three pointings is con- 
sistent with the 1.2 mm continuum distribution shown in 
[Thuma et ai1([200l . The peak of the 1.2 mm continuum 
is at the center of the SW pointing which has the largest 
continuum fraction and smallest beam-scaling exponent 
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Figure 1. Pointings toward M82. Z-Spec's 30" beam (FWHM at 250 GHz) is plotted over the single-dish measurements of the ^-^CO 
J = 2~*-l map from l'i'huma et al.l 1(20001 ). The white star marks the 2.2 ^im peak at oj2000.0 = 9''55™51.9% (5j2000.0 = 69°40'47.14" which 
is the adopted center for the observations. The offsets relative to the center and integration times are given in Table [T] All three pointings 
use a 90" symmetric azimuthal chop throw. 



Table 2 

Continuum Fit Results 



Pointing 


Continuum Fit 


Reduced 




A 


B 


NE 


0.318 ± 0.001 


1.03 ± 0.03 


2.4 


CEN 


0.413 ± 0.001 


1.05 ± 0.02 


3.2 


SW 


0.425 ± 0.001 


0.99 ± 0.02 


3.7 



Note. — The continuum fit columns give the coupling 
fraction (A) and be am scaling exponent {B) of the fit 
described in Section 13.11 and defined in equation [3] The 
reduced values are for the line and continuum fits 
and are calculated based solely on statistical errors and 
do not include errors due to calibration. The fits use 
all data points except the lowest five channels in each 
spectrum and have 132 degrees of freedom. 

of the three pointings. The continuum peak is within the 
CEN pointing and its fraction is almost equal to the SW 
but with a higher index. The NE pointing's continuum 
fraction is the smallest because it is off the peak. 

3.2. Spectral Line Fitting Results 

The channel-to-channel spacing in Z-Spec is 500- 
1300 MHz, corresponding to 700-1200 km , thus the 
instrument does not resolve the line profiles in M82. Nev- 



ertheless, it is possible to fit integrated intensities and 
center frequencies by comparison with the instrumental 
response for each bolometer, carefully measured in the 
laboratory with a long-path Fourier-transform spectrom- 
eter. Each spectral line (indexed by j) is modeled as a 
Gaussian profile of center frequency Vj, amplitude Aj, 
and FHWM 5vj . The sum of these line profiles plus the 
continuum is multiplied by the (normalized) measured 
spectral profile of each channel used in the fitting (in- 
dexed by i) fi{v) and integrated over the range of the 
measured profiles: 180-320 GHz. This process creates a 
model Z-Spec spectrum, which can be compared to the 
observed spectrum. The input line frequencies and am- 
plitudes are then varied to arrive at a least-squares fit, 
using statistical (1/crf) weighting. The frequencies of the 
known lines are fixed in the fitting, modulo a common 
redshift, so that the fit determines the redshift but not 
the line frequencies. In practice for local-Universe galax- 
ies, the ^^CO transition dominates the redshift determi- 
nation. With this method, we obtain accurate centroid 
measurements of line features that did not have an ob- 
vious identification. The fitted line intensities are given 
in Table [3] using the adopted Aw = 250 km s^^ for all 
lines in the three pointings; reducing the FWHM to 100 
km s~^ produces nearly identical integrated intensities 
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Figure 2. Z-Spec 190-307 GHz spectra toward three positions in the nucleus of M82: NE at the top, CEN in the middle, and SW on the 
bottom. The black histogram with error bars are the calibrated measurements and the purple histogram is the spectral fit, including both 
continuum and 20 fitted lines. The fitted continuum is also plotted in dashed dark green and the fitted lines are marked with vertical blue 
and red dashed lines. The dominant CO J = 2—^1 line is shown in gray on the plot, scaled down by a factor of ten; the fit to the CO line, 
also scaled down, is shown in lavender. The CH3CCH J = 16^ 15 transitions and above are shown for reference but are not included in the 
fit. We have identified all three sigma or greater spectral-line features that appear in all three pointings. Nonetheless, there is additional 
structure in the spectra probably due to a multitude of blended, weak lines which Z-spec cannot individually identify. The results from the 
fits are given in in Tables [2l and [3l 



and fit quafity. We comment here on our spectral line 
measurements and how they compare to those found in 
the literature. 

3.2.1. CO and its Isotopologues 

The CO, i^CO, and C^^O J^2^ 1 transitions have 
been exten sively studied and mapped by pr evious ex- 
periments (|Mao et all [20001: iWeifi et al.ll2001[ ): however, 
direct comparison with published intensities is limited 
by beam size mismatch. iWild et al.i (|1992i) quote inte- 



grated main-beam intensities nearly twice what we mea- 
sure in a 13" beam for the CO J = 2 — > 1 transition, 
implying a beamsize (9) scaling of / oc 0~'^' ^ whereas 
the main species intensities given in lMao et al.i (2000) in 
a 22" beam indicate a beam scali ng exponent of —1.2 
to -1.5. The maps obtained by iThuma et all (|2000f ) 
show that the CO emission is more extended than the 
continuum emission which would indicate the exponent 
should be slightly less than unity. The isoto pologues 
13 CO and C^^O have been measured in lMao et all ((20001) 
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Table 3 

M82 Line Fit Results 



Species Sz 


Rest Freq. 


J^uppcr 


Beam Size 


Integrated Line Intensity (K km s ^) 


Transition 


l^^jXlZ ) 




(arcsec) 


NE 


CEN 


SW 


CO J = 2^1 


230.538 


16.6 


32 


423.9 ± 3.0 


460.0 ± 3.2 


452.5 ± 3.3 


i^CO J = 2-i-l 


220.399 


15.9 


34 


29.4 ± 0.7 


31.7 ± 0.7 


30.2 ± 0.7 


ni8n r_o_,i 


21 Q fsfsn 

ij J- c7. ijijyj 


15.8 


34 


8.6 ± 0.7 


7.9 ± 0.7 


7.8 ± 0.7 


riivr T — 9_vi 

V_'l'l t/ ^ r -L 




16.3 


33 


14.3 ± 0.6 


23.3 ± 0.8 


19.9 ± 0.6 


CS J = 4— >-3 


195.954 


23.5 


37 


4.3 ± 0.8 


4.4 ± 0.8 


4.1 ± 0.7 




244.936 


35.3 


31 


^ 1 K ri 01 

< l.D [i.^J 


Z.o it U.O 


/ 1 fl 1 1 

< 1.0 [i.ij 


CS J = 6-i.5 


293.912 


49.4 


26 


< 3.1 [0.8] 


< 3.3 [3.0] 


< 3.1 [1.5] 


riL>W J — o — r Z 


ZD ( .000 


zo. / 


zy 


13.7 ± 0.8 


16.9 ± 0.9 


15.4 ± 0.9 


rlL^lN J — o — r Z 


ZDO.ooD 


zo.u 


zy 


5.1 ± 0.8 


7.0 ± 0.8 


6.9 ± 0.8 


ni>i t/ — o ? z 


971 QJ^I 


26 1 


9R 


2.4 ± 0.6 


2.7 ± 0.7 


3.5 ± 0.6 


C2H J = 3-s-2 


262.251 


25.2 


29 


7.9 ± 0.7 


11.1 ± 0.7 


11.0 ± 0.6 


H2CO (3l3^2l2) 


211.211 


32.1 


35 


4.0 ± 0.6 


2.6 ± 0.6 


< 1.6 [1.5] 


H2CO (3o3-^2o2)* 


218.222 


21.0 


34 


2.3 ± 0.6 


3.6 ± 0.6 


3.1 ± 0.6 


H2CO (4i4^3l3) 


281.527 


45.6 


27 


< 2.5 [1.7] 


< 2.6 [1.6] 


3.3 ± 0.9 


H2CO (404^303)* 


290.623 


34.9 


27 


< 3.1 [1.0] 


3.7 ± 1.0 


3.6 ± 1.0 


H2CO (4i3^3i2) 


300.837 


47.9 


26 


4.0 ± 1.1 


< 3.5 [2.6] 


4.0 ± 1.1 


CH3CCH J = 12^.11 


205.081 


64.0 


36 


5.1 ± 0.6 


4.3 ± 0.6 


4.2 ± 0.6 


CH3CCH J = 13^12 


222.167 


74.6 


33 


2.5 ± 0.7 


3.4 ± 0.7 


2.5 ± 0.7 


CH3CCH J = 14^-13 


239.252 


86.1 


31 


< 1.9 [1.0] 


2.2 ± 0.6 


2.1 ± 0.6 


CH3CCH J = 15^-14 


256.337 


98.4 


30 


< 2.2 [1.0] 


< 2.2 [1.4] 


2.3 ± 0.7 



Note. — Measured line intensities and their estimated uncertainties or the three-cr upper limits for the 
20 identified transitions in the three pointings. For for the upper limits, the value is three times the l-cr 
uncertainty in the fitted intensity and the fitted intensity value is given in square brackets. The adopted 
linewidth for all measurements is 250 km s~^ (sec text); changing this value by factors of two in cither 
direction does not significantly change the fitted line intensities or the quality of the fit. The errors arc 
based on the statistical errors in the spectral data and do not include any uncertainty due to calibration. 
The hyperfine splitting in the CN and C2H transitions is not accessible at Z-Spec's resolution. The H2CO 
lines marked with an * may be blended with other higher energy H2CO transitions. 



and I Wild et al.l (|1992[ ). respectively, and though they in- 
dicate different beam scaling exponents as they did for 
the main species, they both imply that the isotopologue 
distribution is slightly more concentrated than the main 
species. 

3.2.2. CN and C2H 

We detect both the CN J = 2 ^ 1 and C2H J = 3 ^ 2 
transitions in all pointings. The primary energy levels 
of both of these molecules are split by hyperfine inter- 
actions, giving spectra with more structure than is ac- 
cessible using Z-Spec. The intensities quoted in Table 
[3] are obtained by fitting a single Gaussian profile and 
thus represent the total integrated intensity for all tran- 
sitions. These species have be e n det ected, for the first 
time, in M82 by Hen kel et al.l (I1988D and iFuente et all 



(HOOS), but the J = 3 -J> 2 transition of C2H is a new 
detection. Neither of these authors quotes the intensity 
seen in all hyperfine components, making a direct com- 
parison difficult. Our emphasis is on a study of physical 
conditions which, if it included these species, would re- 
quire knowledge of the intensities of the hyperfine com- 
ponents. Therefore, we do not include the CN and C2H 
in this analysis that follows. 

3.2.3. HCO+, HCN, and HNC 

The J = 3 ^- 2 transitions of the HCO+, HCN, and 
HNC molecules are strongly detected in all three point- 
ings. The HCO+ and HCN transitio ns have been previ- 
ously detected bv I Wild et al.l ()1992D . however, the only 
previous detection of HNC in M82 has been the J = 1 — )■ 



line by iHuettemeister et al.l (119951). Compa ring the 12" 
beam measurements from I Wild et alj (| 19921 ) to ours sug- 
gest a beam scaling exponent of roughly —1, consistent 
with the continuum and CO values. We use our mea- 
surements and measurements of other transitions in the 
subsequent analysis. 

3.2.4. CS 

The J = 4 3, 5^4, and 6^5 transitions of CS lie 
in the Z-Spcc band and it's simple ladder of rotational 
transitions make it an ideal candidate for study with our 
instrument. Unfortunately, the lines in our band are not 
very bright in M82 so we can only give upper limits for 
the J = 6 — 5 lines and the J = 5 — >■ 4 lines in the NE and 
SW pointi ngs. The tra nsitions we detect have been seen 
in M82 bvlBavet et al.l (2008) for the CEN pointing and 
IBavet et al.l ()2009f ) for the NE and SW pointings. Longer 
integration times should enable the first detections of 
the J = 6 — > 5 lines, particularly in the CEN pointing 
where the fitted intensity is 2.1a. As with HCO+, HCN, 
and HNC, we combine our measurements with the other 
measured transitions for the analysis that follows. 

3.2.5. H2CO 

We have identified the 3o3— >'2o2, 3i3— 7'2i2, 4i4— >-3i3, 
4o4 — 3o3, and 4i3 — 3i2 transitions of formaldehyde 
(H2CO) which are detected in at least one of the three 
pointings; the latter four of these transitions have not 
been detected previously. Two of these transitions, 
3o3 — > 2o2 and 4o4 — 3o3, may be blended with other 
higher energy H2CO lines and our reported integrated 
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intensity should be interpreted as the sum of the inten- 
sities of all of these blended lines. The feature at 218 
GHz has been studied with the heterodyn e receiver ar- 
ray H ERA on the IRAM 30m telescope by iMiihle et all 
()2007| ). They detected three formaldehyde transitions, 
3o3 2o2, 822 221, and 821 220, near our NE 
and SW pointings along with a possible detection of a 
methanol line in the NE pointing. In addition, they 
demonstrated the power of using formaldehyde lines to 
trace both temperature and density of the molecular gas. 
Higher-resolution followup of the lines we have identified 
would provide a powerful extension to their work. 

3.2.6. CHaCCH 

Methyl acetylene (CH3CCH) has the largest number of 
atoms of any molecule detected outside our galaxy and 
several transit ions have been seen previ ously in the SW 
lobe of M82 (jMauersberger et al.|[l99ll ). Our measure- 
ments in the NE and CEN pointings as well as all the 
measurements of the J=12— J^ll and J=15— >-14 transi- 
tions are new. We suspect that the J = 15 — >■ 14 transition 
may be contaminated with an unidentified feature that 
is causing the fitted flux in the SW pointing to be un- 
realistically high. This spectral feature is right at our 
three-sigma threshold and additional data are needed to 
precisely determine the line identifications. 

4. ANALYSIS 

4.1. Excitation and Radiative Transfer Modeling 

We turn now to a study of the physical conditions in 
the dense g probed with rotational transitions of 
HON, HNC, HCO+, and CS. To make useful inferences 
about the conditions in M82's molecular gas, modeling 
is required. The approach is to adopt basic input pa- 
rameters such as total amount of gas in the beam (col- 
umn density), gas density, temperature, and abundance 
of the species under consideration and then calculate the 
resulting line intensities. A grid of such calculations over 
ranges of input parameters can then provide a framework 
to interpret the observations. We note that the transi- 
tions we study have a range of critical densities ranging 
from 10'* to 10^ cm-^. 

We use the RADEX code (jvan der Tak et al.|[200l ) for 
our excitation and radiative transfer modeling. The pri- 
mary inputs to RADEX are the choice of molecule, the 
kinetic temperature of the molecular gas Tkin, the den- 
sity of molecular hydrogen riHa in cm~'^, and the col- 
umn density of the species iVmoi in cm~^. The radiative 
transfer calculation depends on A^moi/Au, where Au is 
the velocity width of the line. Z-Spec cannot measure 
the linewidth, so /S.v = 250 km s~^ is used through- 
out the radiative transfer modeling to be consistent with 
the linewidth used for the spectral fits. RADEX does 
not assume local thermal equilibrium, but uses an escape 
probability formalism that connects the optical depth to 
the chance an emitted photon escapes the source cloud. 
Several different physical models for this escape proba- 
bility have been derived; we use the expanding spherical 
shell model. However, the results are very insensitive to 
the choice of escape probability. Starting with an initial 
guess for the level population distribution, RADEX com- 
putes the optical depths of all the molecule's transitions, 
from which a new level population distribution can be 



calculated. This process iterates until a self-consistent 
solution is achieved such that the optical depth changes 
by less than a default tolerance from one iteration to the 
ne xt. We use the co llisio nal excitation rates calcul ated 
in lLiaue et all (200^ and lLigue fc Spielfiedell (|2007l ) for 
the CS spe cies and the rates f rom the online database 
outlined in iSchoier et al.l ()2005[ ) for the three remaining 
species. 

4.2. Parameter Likelihood Estimation 

Calculated line intensities are then compared to the in- 
tensity measurements of the species under consideration. 
An additional area filling factor parameter $a must be 
included in the model because the clumps of gas pro- 
ducing the radiation do not in general fill the beam. $a 
scales down the line intensities from RADEX so that they 
can be directly related to the measurements. When $a 
is less than 1, iVmoi represents the column density of an 
individual radiating clump while the product "^A-^moi is 
the beam- averaged column density, <iVmoi>- 

The line-intensity measurements are inherently uncer- 
tain and the physical interpretation using models should 
reflect that. A Bayesian method for calculating likeli- 
hood distributions for various physical quantities of in- 
terest c an be used to address this me asurement unce r- 
tainty [20021: IWard et al.l [20031 [hereafter |Mq1). 

The method constructs the probability distribution of 
obtaining the measurements with their associated errors 
given a set of physical parameters, assuming the mea- 
surements are independent and the errors are Gaussian 
distributed. Using a prior-probability density function 
for the range of physical parameters, the probability dis- 
tribution of the measurements given the physical param- 
eters can be inverted into a likelihood distribution for the 
physical parameters given the measurement results. 

Bayes' Postulate says that the prior probability den- 
sity function should be uniform for all cases in the ab- 
sence of prior knowledge. The prior probability distri- 
butions used for this analysis are assumed to be log- 
arithmically uniform in all model parameters. How- 
ever, this prior probability is used to exclude certain 
non-physical situations relating to large column densi- 
ties and small molecular hydrogen densities. Both con- 
straints require knowledge of the molecular abundance 
ratio, X „i„\ = n mn\ I ?t-Ho . These constraints, described in 
detail by IW08I . limit the total molecular mass contained 
in the telescope beam to less than the dynamical mass of 
th e gala xy ( 2.0 x lO^Mfri, based o n the estimates given 
bv [W03l and iPanuzzo et all (fMol) [hereafter [pTg) and 
limit the column length, equal to the column density di- 
vided by the number density, to less than the length of 
the bright molecular emission on the plane of the sky. 
In addition, models with optical depths larger than 100 
in any transition are excluded because very large optical 
depths are not appropriate for the species under consid- 
eration and RADEX is not accurate when the optical 
depth is this large. 

A critical aspect of the analysis is the scaling of the 
published measurements of various transitions to a com- 
mon beamsize. The limits of beam scaling are (P for a 
source that fills the beam for all measurements, and d~'^ 
for a source that is always smaller than the beam. M82 is 
in an intermediate range with respect to the 25"- 35" Z- 
Spec beam; it is neither fully point-like nor beam-filling. 
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Based on the CO and mm- wave continuum maps, we use 
an intermediate beam scaling of 0^^, appropriate for the 
distribution which is to first order extended along the 
major-axis, but unresolved along the minor axis. 

4.3. Multi-Species Model 

Four of the species detected in this survey have the re- 
quired radiative and coUisional rate data available in an 
online database dSchoier e^ al. 2005) or in the literature 
(|Liaue et al.l 120061 : iLiaue fc Spielfiedeil 120071 ): HC O+, 
HCN, HNC, and CS. The method described in lWOl and 
outlined in the previous section was developed for analyz- 
ing 11 transitions of CO and ^'^CO and it can be applied 
to each of the four molecules above individually. How- 
ever, for three of the four selected molecules, the number 
of detected transitions (see Table |4]) is less than the four 
primary parameters of the radiative transfer model (Tkin, 
'^Ha: ^moi, and $a)- CS and C^'^S have had several tran- 
sitions measured but with relatively low signal-to-noise 
which would lead to very broad constraints from the like- 
lihood analysis. 

Observations of these four species in star-forming re- 
gions within our Galaxy indicate good correspondence 
with both the spatial distributions and t he line pro- 
files (|Brand et al.|[200lt iNikolic et al.|[200l . High spa- 
tial resolution maps of t he Galactic circumnuclear disk 
(jChristopher et all I2005D and of the starburst galaxy 
NGC 253 (|Knudsen et all [20071) in HCO+ and HCN in- 
dicate strong similarity, in general, between the emission 
of these two molecules. 

We thus construct a model in which all four species are 
characterized by a common kinetic temperature, molecu- 
lar hydrogen density, molecular hydrogen column density 
{Nh^), and filling factor. Each species is modeled with a 
individual abundance (X,„oi). Extraction of the absolute 
abundances relative to H2 is not possible with RADEX, 
but it can constrain the relative abundances of the vari- 
ous species. In the following analysis, CS is chosen as the 
primary species and the model is parameterized by Tkin, 

J^Ha: -^CS, '^'A, -'^^HCO+Z-'^CS, -'^HCn/-'^CS, -^HNc/-'^CS, 

and Xc34s/Xcs- 

The result of the likelihood analysis is a likelihood ma- 
trix with each point in the matrix characterized by a 
particular value of the four primary species parameters 
and three secondary species abundance ratios whi ch are 
used to parameterize the model. As described in IW031 
likelihood distributions for a single parameter can be ob- 
tained by integrating the likelihood matrix along all the 
other dimensions. These distributions, seen in Figures |3] 
and m can be used to calculate both median values and 
confidence ranges for the seven parameters in our model 
(see Table [S]). It is also possible to calculate likelihood 
distributions for parameters which are functions of model 
parameters, such as the gas pressure, P — x Tkin, and 
beam-averaged column density, <A^cs>= '^axA^cs- The 
beam-averaged column density can be used to calculate 
the total molecular mass in the beam by 



1.5 X m^f, X 



<A^mol> 



(4) 



where m/fa is the mass of a hydrogen molecule, -Dbcam is 
the linear diameter of the beam in cm and the factor of 
1.5 accounts for the additional mass of He and dust in the 
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Figure 3. Likelihood distributions for the four primary species 
parameters, T]jin, nn^, Nqs, and <3?Ai for the five-species radiative 
transfer likelihood analysis. All of the distributions have well de- 
fined peaks, indicating the parameters are well constrained by the 
measurements within the range of the RADEX simulation. The 
distributions for the three pointings are plotted in dash-dotted, 
solid and dashed lines for the NE, CEN, and SW pointings, re- 
spectively. The distributions are normalized to have unit integral 
when integrated over the base- 10 logarithm of the parameter value. 
The jaggedness seen in the Tkin distributions is caused by two fac- 
tors: first, the spike seen in the distributions for all three pointings 
at 300 K is cau sed by p atching together the two sets of colisional 
rates of CS from'Lique et al. ( 2006 ) and jLigue fc SpielfiedeL (2007). 
The former paper calculated the rates for 31 rotational levels for 
temperatures up to 300 K while the later produced rates for 38 
rotational and 3 vibrational levels in a higher temperature range. 
The remaining jaggedness in the CEN pointing Tkin distribution 
probably comes from a numerical problem in RADEX where, for 
certain physical conditions, it fails to converge on a stable solution 
for the optical depth in the CS J = 1 -> line. The NE and SW 
pointings' distributions do not show this effect because they do not 
have a CS J = l— ^-O measurement. 



molecular clouds. Distributions for these parameters and 
their associated two-dimensional distributions are shown 
in Figure O 

The emitting regions containing these molecules are 
likely to be isolated cores of high-density gas. In that 
limit, we can estimate the velocity gradient inside the 
cores with a simple geometrical argument. The charac- 
teristic size of an individual core can be estimated by 
•S'corc « (-^Vcs/-^cs)/"-H2, which is simply the length de- 
fined by the ratio of the molecular hydrogen column and 
volume densities. This value can be used to estimate the 
number of cores in the beam by computing the ratio of 
the area of the emitting region in the beam over the size 
of a single core, iVcoro « (*A x -Dboam)/'5'coro- A reason- 
able approximation for isolated cores where $a <C 1 is 
that the observed total line width Av is split up equally 
among the individual cores. That implies that the veloc- 
ity gradient in a single core is 



dv 
dr 



Av/N„ 



Av 



CS 



$A?^H2 



(5) 
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Table 4 

High-Dipole-Moment Species Measurements 



Transition 


Rest Freq. 


-^uppcr 


Obs. Beam NE Flux 


CEN Flux 


SW Flux 


Refs 






(GHz) 


(K) 


(arcscc) 


(Iv km s ) 


(K km s-i) 


(K km s-i) 
















cs 










J = 




48.991 


2.4 


36 




16.2 ± 1.1 




2 




J = 


2-i«l 


97.981 


7.1 


25.1 


9.4 ± 0.2 


13.3 ± 0.3 


8.9 ± 0.2 


3. 


4 


J = 


3-i-2 


146.969 


14.1 


16.7 


8.9 ± 0.1 


11.2 ± 0.3 


7.6 ± 1.7 


3, 


4, 5 


J = 


4->3 


195.954 


23.5 


37.1 


4.3 ± 0.8 


4.4 ± 0.8 


4.1 ± 0.8 


1 




J = 


5-S-4 


244.936 


35.3 


30.8 


< 0.5 


2.3 ± 0.5 


< 0.5 


1 




J = 


6-i>5 


zyo.y iz 


ACi A 


26.3 


< 1.0 


/ 1 1 
< i.l 


< l.U 


1 
























J = 


3->2 


144.617 


13.9 


17 


~0.5 ± 50% 


0.6 ± 0.1 


~0.4 ± 50% 


6 




J = 


4-5.3 


192.818 


23.1 


37.1 


< 0.8 


< 0.8 


< 0.8 


1 




J = 


5-5-4 


241.016 


34.7 


30.8 


< 0.5 


< 0.5 


< 0.5 


1 




J = 


6-5.5 


289.209 


48.6 


26.3 


< 1.0 


< 1.1 


< 1.0 


1 














HCO+ 










J = 


1-^0 


89.189 


4.3 


23 


35.0 ± 2.0 


38.3 ± 2.0 


37.2 ± 2.0 


7 




J = 


3^2 


267.558 


25.7 


28.6 


13.7 ± 0.8 


16.9 ± 0.9 


15.4 ± 0.9 


1 




J = 


4-5.3 


356.734 


42.8 


14 


23.6 ± 1.7 


22.2 ± 1.7 


22.6 ± 1.9 


8 














HON 










J = 


1-^0 


88.632 


4.3 


23 


21.8 ± 2.0 


18.5 ± 2.0 


23.7 ± 2.0 


7 




J = 


3^2 


265.886 


25.5 


28.8 


5.1 ± 0.8 


7.0 ± 0.8 


6.9 ± 0.8 


1 




J = 


4->3 


354.505 


42.5 


14 


5.6 ± 0.6 


9.0 ± 0.7 


6.1 ± 0.4 


8 














HNC 










J = 


l->0 


90.664 


4.4 


25 


10.7 ± 2.0 


13.4 ± 0.8 


12.3 ± 3.0 


9 




J = 


3-5.2 


271.981 


26.1 


28.2 


2.4 ± 0.6 


2.7 ± 0.7 


3.5 ± 0.6 


1 





References. — (1) This Wor k; (2) Paelionc ct al. 1995; (3) Bayct ct al. 2008; (4) Bayct ct al. 2009' ; (5) 
I Mauersberger &€ Henkel||198g|: (6) [Martm ct al. 2009; (7) Nguycn-Q-Ricu ct al. 198a ; (S) .Scaquist &; Frayed 
I2OOOI : {"Ql lHuettemeister et al.lll995[ 

Note. — Fluxes and upper limits used in our multi-species radiative transfer likelihood analysis. In 
addition to the given statistical error, a 10% calibration error is added in quadrature to each measurements' 
uncertainty. An additional 10% error is added to the lines measured in beams smaller than 18'' and to the CS 
J — 1 — > line. The CS -7—1—^0 line has only been measured for the CEN pointing and thus is not included 
for the likelihood analysis for the other pointings. The C^'^S J — 3 —¥ 2 has also only been measured in the 
CEN pointing; we estimate the flux in the other pointings based on the J — 3^2 line ratio between CS and 
C'^^S in the CEN pointing and apply a 50% error to these estimates for the likelihood analysis. The C^^S 
upper limits are based on the uncertainty in the spectral fit for the main-species lines; the actual integrated 
intensity in the C^'^S lines is well below the detection threshold of our data. 



Table 5 

Mult i- Species Modeling Results 



Quantity 


NE Pointing 


CEN Pointing 


SW Pointing 




Median Range 


Median Range 


Median Range 



Primary Species Parameters 



Tkin (K) 


160 


58 - 


470 


130 


56 - 320 


130 


44 ■ 


- 410 


"•H2 (cm ^) 


104.3 


104 _ 


104-7 


104-2 


104.0 _ 104.5 


104-4 


104-1 


- 104-8 


A^CS (cm-2) 


1015.7 


1015.5 _ 


1015.9 


IOI6.O 


IOI5.9 _ 10I6.2 


IOI5.7 


1015.5 


- 1015-9 


<S>A 


10-21 


10-2.3 _ 


10-1.9 


10-2-1 


10-2-^ - 10-2-0 


10-2-1 


10-2.3 


- 10-2-0 








Second 


ary Species Relative A bund 


ances 








lOO.i 


lOO O - 


IOO.2 


IOO.O6 


10-0.02 _ 100.13 


IOO.2 


lOO.i 


- IOO-3 


^hcn/-''^cs 


100-4 


100-3 _ 


IOO.5 


IOO.3 


IOO.2 _ 10O.4 


IOO.5 


IOO.4 


- IOO-6 


^HNC/'''^CS 


IqO.O 


10-0.2 


- lOO-i 


lOO.o 


10-0.2 _ loO.l 


IOO.2 


IqO.O 


- IOO-3 


Xc34s/Xcs 


10-1-^ 


10-2-7 - 


10-1.3 


10-1.6 


10-1-7 - 10-1-4 


10-1.8 


10-2-7 


- 10-1-4 










Projected Parameters 








Pressure (K cm~^) 


106.6 


106-3 _ 


106.8 


106.4 


106.2 _ 106.6 


106.5 


106.3 


- 106-« 


<Ncs> (cm- 2) 


1013.6 


1013.5 _ 


1013.8 


IOI3.9 


IOI3.8 _ 1014.0 


IOI3.6 


IOI3.4 


- 1013-7 


dv/dr (km pc" 


1.9 


0.6 - 


4.4 


4.5 


3.0 - 7.2 


1.6 


0.5 


- 4.1 


Total Gas Mass in Beam (Mq) 


10T.9 


107.8 _ 


108.1 


IOS.2 


108.1 _ 108.3 


107.9 


107-8 


- 108-0 



Note. — Results obtained from the multiple species radiative transfer modeling of the lines of CS, HCO"'', HCN, HNC and 
C^'^S. The median and 68% (lo") confidence ranges are obtained from the likelihood distributions shown in Figures fsl — fsl 
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Figure 4. Likelihood distributions for the abundances of the sec- 
ondary species, HCO+, HCN, HNC, and C^*S in blue, red, green, 
and dashed black, respectively, for the NE (top), CEN (middle), 
and SW (bottom) pointings. The thick vertical black line indicates 
the solar abundance of '^''S/^^S = 0.043 that is assumed for calcu- 
lating all other likelihood distributions. The bottom axis is labeled 
with the abundance relative to CS, v^hich is the parameterization 
used in the model. The top axis shows the abundance relative to 
H2 using the adopted value of Xqs = 3 x 10~^. 



The distributions for the core velocity gradient are shown 
in Figure [S] along with two-dimensional projections of the 
three-dimensional distributions used for the calculation. 

4.4. Molecular Abundances 

The likelihood distributions of the HCO+, HCN, HNC, 
and C'^'^S abundances relative to that of CS are shown in 
Figure m These represent the first statistically rigorous 
measurements of molecular abundances in M82. Aver- 
age values for the entire starburst nucleus are presented 
in Table [S] As discussed in Section the abundance 
of CS is used in the likelihood calculation to apply cer- 
tain physical limits on the parameter space of the radia- 
tive transfer grid. We adopt the CS abundance Xgs o f 
3 X 10"^ calculated in IMauersberger fc Henkel (| 19890 . 
Their calculation is based on assuming optically thin 
CS emission and the CO intensity to H2 column den- 
sity conversion factor that is observed in the Milky Way. 
Changing the CS abundance by half an order of mag- 



Figure 5. The plots in the left column show two-dimensional like- 
lihood distributions for three pairs of primary species parameters 
while the right column plots show distributions for parameters that 
are projected from the corresponding 2D distributions. The three 
pointings, NE, CEN, and SW, are indicated by red, black, and blue 
contour lines in the 2D distributions and with dash-dotted, solid, 
and dashed lines in the projected ID distributions. The contour 
lines represent 10%, 20%, 40%, and 80% of the peak of the 2D like- 
lihood distributions. The top-left plot has dashed lines of constant 
pressure, labeled in units of logjQ K cm~'^; the Tj^in versus rin^ 
distributions are used to create the pressure distributions shown 
in the top-right plot. Dashed lines of constant beam-averaged col- 
umn density in the left-middle <I>a versus Nqs plot show how the 
<^CS> distributions in the right-middle plot are obtained. The 
alternative axis in the beam-averaged column density plot shows 
the molecular mass in the beam assuming Xcs = 3 X lO"'*. The 
bottom-left plot shows the N^q versus rin^ distributions; the diag- 
onal dash-dot-dot-dot lines indicate constant velocity gradient in 
units of km s~^ pc~^ at the most-likely value of <I>a- The velocity 
gradient distributions shown in the bottom-right panel are com- 
puted from the three-dimensional distributions of Nqq, n^^, and 



nitude either up or down does not affect the likelihood 
distributions for nearly all of the parameters; the molec- 
ular hydrogen density and kinetic temperature distribu- 
tions shift down and up, respectively, with increasing CS 
abundance such that the distribution of gas pressure is 
relatively unchanged. The agreement between the mea- 
sured and most-likely model's integrated line intensities, 
shown graphically in Figure |6] for Xqs = 3 x 10~^, is not 
significantly impacted by changing Xqs- Furthermore, 
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Table 6 

Molecular Abundance Ratios for M82 



Species 


CS 


HCO+ 


HCN HNC 


C34S 


loSlo(^mol) 


-8.5 


-8.4 


-8.1 -8.5 


-9.9 



Note. — Abundance ratios (with respect to H2) de- 
rived from the multiple species modeling results for the 
abundance of HCO+, HCN, HNC, and C^*S relative 
to CS and the abundance for CS relative to H2 from 
IMauersberger &: Henkell il989l) . The uncertainty in the 
four modeled abundances is ±0.1 and is less than ±0.5 
for CS in logarithmic units. 



while the hkehhood distributions for the C'^^S abundance 
show shght differences between the three pointings, we 
have httle reason to suspect that this value would be 
much different from the solar isotopic abundance ratio 
of '^''^S/^^S = 0.043. Therefore, we impose this isotopic 
ratio for the other likelihood calculations. 

Of course, changing the assumed CS abundance im- 
pacts the conversion from the relative abundance ra- 
tios to standard abundances relative to H2. Also, while 
the beam-averaged column density distribution does not 
change when using different values for the CS abun- 
dance, the conversion from <Ncs> to total mass is in- 
versely proportional to the abundance (see equation 2]). 
Of the thre e values of CS abundance we h ave used, the 
value from IMauersberger fc Henkell (|1989f ) seems to be 
the most reasonable; lowering it drives the total molec- 
ular mass too high and raising it increases the tempera- 
tures to extreme values. 

While our likelihood analysis is an independent 
confirmation of the CS a bundance estimated by 
IMauersberger fc Henkell (|1989[ ). it should be noted that 
there is significant debate about the CO intensity (/co) 
to molecular hydrogen column density (iViJa) conversion 
factor. Acq = Nn^/Ico^ used in their calculation of 
CS abundance. Their value was 2.2 - 2.5 times larger 
than more recent measurements of Xgn in the Milky 
Way (1.6 X lO^^ cm-^ (K km s- ^)-\ iHuiiter et al.|[T997l 
1.8 X 10^0 cm-2 (K km s~M~ MDa^e et al.l[200TI ). Fur- 
thermore, a detailed analysis of an interferometric CO 
J = 1 — map of the nucleus of M82 showed that 
Xco is roughly 2- 10 times less tha n what is measured 
in the Milky Way (|WeiB et al.ll200lD . Quantitative mod- 
eli ng of ensembles of unres olved giant molecular clouds 
bv lMalonev fc Blackl (|1988D showed that Xco IS propor- 

tional to T^^^ ; this proportionality was confirmed 
in M82 by iWeifi et all (|200H ). Studies of several galax- 
ies have also shown th at Xgn also depends on metalic- 
ity (Wilson 1995; Bos elh et all [2002^ and CO intensity 
iNakai & Kunol I1995D . The deviations from the stan- 
dard M ilky Way value can be significant; .Casasola et alj 
()2007[) measured the conversion factor in several giant 
molecular associations in M81, the i ntera ction partner 
of M82 and NGC 3077 (lYun et al.l[T99l . and found 
the galactic average of Xqo to be 20 times larger than 
the Milky way value. In spite of these developments 
in the understanding o f Xgo since the publication of 
IMauersberger fc Henkell ()1989[) . our results indicate that 
their estimate for the CS abundance is correct to within 
a half an order of magnitude. 
Adopting the value Xqs = 3 x 10^^ allows for compar- 



ison of our abundance measurements to those found in 
the literature. For the most part, M82 is assumed to have 
abundances similar to those found in regions of high-mass 
star formation in our Galaxy such as Sgr B2 or Orion. 
IWild et al.l (|1992D quote X^co+ = 1 x lO"** and Xhcn = 
2 X 10~^ which have also been used by other authors (e.g., 
ISeaauist fc Fraveri [2000l ). These are somewhat higher 
than what we measure tho ugh we agree that HCN is 
twice as abundant as HCO+. IHuettemeister et all H1995[ ) 
use their measurements of HCN J = 1 — >■ and the CO 
intensity to H2 column density conversion factor to calcu- 
late Xhnc = 1-4 X 10~^°, assuming the HCN emission is 
optically thin. Our measurements of the HNC J = 3 — )■ 2 
contradict the optically-thin assumption and our models 
point to an abundance 1.5 orders of magnitude larger, 
roughly equal to the CS abundance. 

4.5. Mass of Dense Gas 

The total mass of molecular gas in the nuclear region 
can be estimated from the total mass in the beam mea- 
sured for each of the three pointings by adopting a par- 
ticular geometry for the central region. If the area of the 
nucleus is and the area of overlap between our 30" 
beam and this nuclear area is Aqi, then the total mass in 
the nucleus can be calculated using 

Mtotal = A^CEN + MnE ( ^"2^ ^°' ) 

where Mcen, A%e, and Mgw are the masses measured 
in the three pointings. This equation uses the mass sur- 
face density implied by A/ne and Afgw to account for the 
mass not covered on the left and right sides of the CEN 
pointing. If we adopt a rectangular nuclear geometry 50" 
X 15", as indicated by the CO J = 2— >1 interferometer 
maps m IWeifi et al.l (|200l . then the area scaling factor 
{A,,-Aoi)/2Aoi = 0.37 and A/total = 2.2±0.5x 10*^ A^^g in 
the nuclear starburst. This value is comparable to the to- 
tal mass traced in CO from both multi-line studies (W03l : 
2.OxlO*Af0 when corrected to d = 3.86 Mpc) and via 
interf erometric low-J CO and ^''CO imaging (jWeifi et al.l 
I2OOI . 

Since we measure a gas mass with a suite of tran- 
sitions, we can provide a calibration of the HCN X- 
factor, Ahcn, which converts HCN line luminosity (in 
temperature units) to mass of dense gas. Assuming that 
the source couples to the various beam sizes as 6^, we 
find A'hcn.j=i of 10-15, in agreement with the canon- 
ical value of 10 derive d from virial considerations by 
IGao fc Solor^ (pOOlal lbh. Since HCN is sub-thermally 
excited, using the J = 3 — 2 transition Ajjcn is much 
higher with values of 36-65. 

5. DISCUSSION 

5.1. Physical Conditions and Relationship to 
CO-Traced and Atomic Gas 

Our likelihoods suggest temperatures between 50- 
500 K, bro adly c onsistent with the warm components 
modeled by 'WOS, up to J = 6 — )■ 5, an d somewhat lower 
than the 400-800 K derived by iPld in considering all 
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Figure 6. Measured and modeled spectral energy distributions for the five molecules of the multi-species radiative transfer likelihood 
analysis for the three pointings, NE (left), CEN (center), and SW(right). The measurements of CS, HCO+, HCN, HNC, and C^-'S, from 
Table |4] scaled to a common 30" beam, are marked with red stars, brown squares, purple triangles, teal triangles and open blue stars, 
respectively. The error bars include both measurement and adopted calibration error, if any. Three-sigma upper limits are indicated with 
downward arrows using the same color scheme as the measurements. The most-likely model from the likelihood analysis is indicated by 
the six-pointed stars connected by lines, in colors to match the measurements. 



of the transitions up to J = 13 — 12 as measured with 
Herschel SPIRE. Our results are also consistent with the 
temperature of 200 K inferre d from the formald ehyde 
measurements and analysis of iMiihle et al.l (|2007| ). The 
relatively low precision with which we measure the tem- 
perature is not surprising since the most likely tempera- 
tures are generally higher than the upper level energies of 
the transitions we are studying (e.g., HCN J — 3: T = 25 
K) . We do note that our results do not support the pres- 
ence of substantial amounts of gas at temperatures below 
'-^SO K unless the density is so low that the transitions in 
our analysis would not be excited. This means that the 
drop in line intensity with J (in temperature units) is 
due to sub-thermal excitation of the levels above J = 1, 
and the run of line intensity with J should thus provide 
a reliable density measurement. Indeed, our derived me- 
dian densities are 1.5-3 xlO'*cm~'^, sub-critical for the 
transitions above 200 GHz. 

Our derived densities are larger than those der ived 
with CO studies. IWOSi find 600-6000 cm'^ and IpTC 
find 1000-13000, though we do note in some cases IWOS 
find densities poorly constrained on the high-density end. 
The product of temperature and density is the thermal 
pressure, and we find values of 1-4 x 10^ K cm~'^, an order 
of magnitude higher than the pressure inferred bv IW03I 
(0.5- 4 X 10^ Kcm~'^), but comparable to that derived by 
iPld using all of the CO transitions. Our high densities 
may reflect the fact that we are probing preferentially 
high-density cores, in approximate pressure equilibrium 
with the larger, more diffuse envelopes which produce the 
bulk of the CO. This should not be surprising since HCN, 
HNC, and CS are generally found in UV-shielded cores as 
their dissociation energies are less than that of CO, and 
they don't have generally achieve sufficient column densi- 
ties to self-shield. For instance, the photo- dissociation re- 
gion (FDR) chemical models of lFuente et al.. (.2008 ') show 
that HCN exists primarily within Ay > 5. 

We compare our results wi th the studies of the 
photo-dissociated atomic gas. iKaufman et al.l (|1999f ) 



and iColbert et al.' (199 91) have applied a PDR model 
(updated from lWolfire etaLlfTOOOl: iTielens fc HollenbachI 
119850 to [C II] and [O I] fine-structure-line measure- 
ments from the Kuiper Airborne Observatory (KAO) 
and ISO Long- wavelength Spectrometer (LWS), respec- 
tively. PDR conditions are parametrized in terms of the 
density and UV field strength Go. Estimates for the 
M82 central starburst range from: 1) rtH2=10* cm^"^, 
Go=103-5 (IKaufman et al.l [l99l assuming the [C II] 
emission is uniformly distributed over its 55" beam, so 
only a small fraction arises in the [O I]-emitting re gion) , 
2) nH2=10^-^cm-3, Go = 102-5 (|Kaufman et al.l [l999L 
assuming that all the large-beam [C II] emission arises 
in the same r egion as the [O I]), a nd 3) nH2 = 10'^ '^ cm""^, 
Go = 102-« (IColbert et al.l[l999l using large-beam ISO 
fiuxes, but removing a [C II] contribution from the ion- 
ized gas). These estimates form a locus in the nnj, 
Go plane, with the upper end of the density range be- 
coming consistent with our density likelihood. The UV- 
illuminated surfaces of clouds might be expected to have 
somewhat lower density than the UV-shielded cores, and 
we note that with the modeled surface temperatures of 
^300 K, the PDR thermal pressures range from 10^ '^ to 
10^'5, broadly consistent with the values derived from the 
CO as well as our analysis. 

However, the PDR models which fit the atomic line 
fluxes cannot explain the strength of the mid- J CO tran- 
sitions. For the range of PDR conditions inferred from 
the atomic lines, the modeled CO J = 7 — 6 to [C II] 
intensity ratio is at most ~ 1.7 x 10^"^ (at n^lO* cm^'^). 
The CO spectrum indeed peaks (in energy unit s}^ at 
J=7— ^6 per the Herschel SPIRE measurements (|P1Q) . 
This transition carries a fraction 1.4-2.3% of the [C II], 
depending on whether the [C II] is resolved or unre- 
solved in its 55" beam, relative to the 43" CO analy- 
sis region. Thus the CO emission in M82 exceeds the 

^ see also http://dustem.astro.iimd.edu/pdrt/index.html 
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PDR predictions for the measured densities by an order 
of magnitude. Moreover, this is a lower hmit — if the [C 
II] includes a contrib ution from ionized gas (as assumed 
bv lColbert et all [1999) . then the inferred PDR line ra- 
tio is even larger, and less consistent with the model. 
Higher-density PDR models n > lO'*'^ cm"^ can bring 
the mid-J CO to [C II] into agreement with the obser- 
vations, but then the [O I] transition is over-predicted 
relative to the observations by a factors of 3-10. More- 
over, such high-densities for the PDR would be inconsis- 
tent with our measured density likelihood which excludes 
n > lO-^-^cm-^. 

Empirically, the suite of [C II], [O I], and mid- J CO 
emission in M82 does not compare with the PDRs as- 
sociated with Galactic star-formation regions. In M17, 
the bright mid-J CO and its widespread distribution 
with [C II] over several parsecs as viewed in t he ed ge- 
on PDR led lStutzki et all (|1988| ): iHarris et all (|1987D to 
a clumpy PDR interpretation. However, their data indi- 
cate that the CO J = 7 — >■ 6 in this source only amounts 
to some 0.5-1 X 10"'^ relative to the [C II], much less than 
in M82. The Orion PDR is more extreme, wit h CO 
J = 7 ^ 6 some 8 % of the [C II] (Schmidt-Burgl^lMS; 
iStacey et al.lll993l ). but in this case the [O I] 63 /im line 
dominates the PDR lin e emission with 11 x the [C II] (see 
iHerrmann et al.lll997D . so the CO J = 7 ^ 6 is ~0.8% 
of the atomic line emission. While the partitioning of 
energy between the atomic and warm molecular compo- 
nents in the Orion bar is thus similar to what we observe 
in M82, we stress again that the line ratio are not a good 
match, particularly the [O I] to [C II], which is only 1- 
1.5 in M82. The bright CO and [O I] in Orion is consis- 
tent with a clumpy PDR model wh ich includes clumps 
with d ensity as high as 10^ cm~'^ d Burton e t ahl 119901 : 
iKoest er et al. 1994; Mei xner fc Tiele ns 1993). Such den- 
sities cannot be commonplace in M82 based on line ra- 
tios in the atomic gas, the CO analyses, and our analysis 
of the high-dipole-moment species which indicate typical 
densities less than lO^'^cm"'^. 

5.2. Heating of the Gas 

The poor match to the Galactic PDRs and the inabil- 
ity of the PDR models to explain the powerful mid-J 
CO emission suggests that non-UV heating sources may 
be dominating the energetics of the molecular material in 
this starburst nucleus. X-rays can be a powerful source of 
energy input to the gas, produci ng lumin ous X-ray Disso- 
ciation Regions (XDRs. Malo n ev et al.l[T996,) . However, 
the hard X-ray luminosity of M82 is only 1.1 x lO^L© 
([Strickland fc HeckmaiJ[2007l ). completely inadequate to 
power the observed CO emission. Moreover, multiple 
chemical / excitation studies show that the line emission 
from M82 is not consistent with X-r ays being a dominant 
heatin g ter m. The XDR models o f iMeiierink fc SpaansI 
([20051 ) and iMeijerink e t~aD (120071) predict more [O I] 
than is observed ([Colbert ct al. 19 9^ at o ur me asured 
densities around 10 cm" iFuente et al.l ([2008') stud- 
ied HOC+ and HCO+ in detail and showed that the 
HOC+ emission is not spatially correlated with X-ray 
emission, and the line ratios among these species as well 
as C0+ and CN are well-matched with PDR models with 
n > 10* cm~'^ and Go = 10*, similar to the values derived 
in analysis of the atomic lines. They did not consider the 



total energy budget of the molecular gas, and did not dis- 
cuss heatin g mechanisms. 

Simiarlv. iLoenen et al.l ([2008D have considered XDR 
models and PDR models with and without extra bulk 
molecular heating and have predicted line ratios among 
the HON, HCO+, and HNC J = 1 ^ transitions. Ac- 
cording to their models, the measurements for M82 (Ta- 
blelU are inconsistent with XDRs, as XDRs are predicted 
to have HNC comparable to or brighter than HCN. The 
line ratios are best fit with PDRs with n ^ 10*'^cm~'^ 
(their lowest density considered), but with additional 
bulk heating on the order of 3 x 10~^® ergs~^cm~'^, 
^ 1Lq/Mq input into the molecular material (they as- 
cribe this to mechanical heating). 

We thus look for other p otential b ulk heating 
sources other than X-rays. iSuchk ov et all ^93) pro- 
posed that the elevated level of cosmic rays due to 
the supernova rate in M82 will heat the gas, and 
their derived cosmic ray enhancement factor relative 
to the Galaxy of ~500 has been confirmed by re- 
cent measurements of high-ene r gy ga mr na-rays in M82 
(VERITA S Collaboration eFall (20091) . Bradford et all 
( 2003) showed that the energy input from a similar cos- 
mic ray density is indeed well-matched to the total molec- 
ular gas cooling in NGC 253 as extrapolated from the CO 
transitions up t o J = 7— »-6, where the CO em ission ap- 
pears to peak, ([Hailev-Dunsheath et al.l[2008[ ). Because 
the CO excitation and total cooling in M82 are similar 
(from WOl, [Pl3 fluxes), we conclude that cosmic rays 
are a plausible means of heating the gas. 

Another important heating term for bulk molec- 
ular material is the dissipation of turbulence 
felgaronc & Pugct 1 995tlMac Lowll999l : [Fan fc Padoanl 
[20091: Bradford et al]]2005l) . The average heating per 
gram is given by dimensional analysis 0.5 x ^/Sa^/L, 
where L is the driving scale for the turbulence and ay 
is the (1-D) RMS ve locity spread on this size scale. 
[Pan fc Padoanl (|2009[ ) conclude that turbulent dissi- 
pation with ay ~ 2.5kms~^ on 0.1-1 pc scales is the 
dominant heating term on average in Galactic clouds 
(exceeding cosmic ray heating by factors of 3-4), and 
that it produces temperatures of 13-36 K in Galactic 
clouds with densities of 10*-10^ cm~^. For gas densities 
above 10^ cm^'^, gas-grain cooling becomes dominant, 
and limits the increase in the gas temperature, as the 
dust energy balance is largely independent of the gas 
temperature. A lower bound to the total cooling per 
mas s in M82 is given by the total C O lum inosity in 
IPIQ I ratioed to the mass derived by [W03l using the 
mid-J lines, ~0.1 Lq/Mq. Achieving this via turbulent 
dissipation requires ~5-20kms~^ on the same 

0.1-1 pc scales. The velocity gradient implied by our 
radiative transfer analysis of 4-10 kms^^pc"^ is a bit 
lower than this would suggest, but it becomes consistent 
if the turbulence is distributed on few-pc scales (e.g., 
ay ^ 20kms^^, L ^5 pc), as might be expected if it 
is produced by winds and supernova shocks from young 
star clusters. 

5.3. Comparison with Expanding Shell Starburst 
Models 

Our findings are largely consistent wit h th e evolv ing 
starburst model of .Yao et al.. (,2006. ) and[^ ([20091) . in 
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which the gas in the nucleus of M82 is a superposition of 
expanding spherical bubbles around stellar clusters. The 
bubble interiors are ionized gas ([H II] regions), and the 
shells are swept-up molecular gas, so the inside edges are 
PDR fronts. The PDR surfaces are responsible for most 
of the atomi c line emission, and would agree with the 
iFuente et all (j2OO80 resuhs. The PDR shell also gener- 
ates most of the excited CO emission, though their model 
includes mechanical energy input from the shock due to 
the expansion, which may be a substantial term in the 
heating of the warm molecular gas. 

Model intensities for the high-dipole moment species 
are not presented, preventing a detailed comparison with 
our data, but (not surprisingly) the total mass and phys- 
ical conditions in their modeled shells are similar to what 
we find with our likelihood analysis. Their estimated gas 
mass of 2 X 10^ Mq from CO in the central 1 kpc along 
the major axis (nearly identical to our modeled 50" x 15" 
region) is comparable to our measured 2.2 x 10^ Mq, and 
the shell density at the putative 3-10 Myr age is modeled 
to be l-3xl0 "''cm~'^, similar to what we find. However, 
as lYaol (|2009f ) notes, there are some inconsistencies in the 
model. While all of the observed molecular and atomic 
line emission is reproduced, the stellar luminosity which 
is required is only ^5% of the observed far-IR luminos- 
ity in the same region, potentially the result of assuming 
zero pressure for the ambient ISM which results in more 
mass swept up in the modeled shells than is physical. 

5.4. Warm Star-Forming Gas 

Regardless of the details, the I Yaol ()2009() model is rep- 
resentative of the likely physical situation: a new stellar 
cluster subjects the surrounding molecular gas to both 
UV photon and mechanical energy input which heats 
and compresses it, at least in the first 10 Myr after the 
starburst. The result is a molecular ISM that is demon- 
strably warmer than the Galac tic cloud cores. Doe s this 
mean that the SF is quenched? iFuente et al.l (|2008[ ) com- 
pare line ratios of ions (HOC"'", CO"*" to HCN) with their 
PDR chemistry model to estimate the total depth of the 
PDR (ions except HCO+ quickly become less abundant 
with increasing ^v)- Fitting the line ratios to two com- 
ponents, they find some ^^87% of the molecular gas is 
in small clouds with Ay ~ 5 (but large enough to house 
HCN, HNC) with only ~ 13% in clouds with Ay ~ 50 
and conclude that in general the molecular gas is highly 
fragmented with clouds too sm all to form massive stars. 
iForster Schreiber et al] (|2003l ) found that the SF in M82 
has occurred in two bursts, one in the center some 10 Myr 
ago and one in a circumnuclear ring ~5 Myr ago, and 
that each burst was self-quenching with a timescale of 
a few Myr due to mechanical energy input into the gas. 
The total mass of stars formed through both episodes 
is modeled at 2-5xlO^M0, depending on the low- mass 
part of the IMF (and cannot be more than ^ 6 x 10^ Mq, 
the total measured stellar mass in the system) . The stel- 
lar mass formed in the last 10 Myr is thus comparable to 
or at most double the amount of dense molecular gas re- 
maining, so unless the eventual star formation efficiency 
is limited to 30-50%, one may ask if the gas can be the 
raw material for another round of SF. 

If the material we trace is indeed forming stars, then 
the warm molecular medium is likely to impact the stel- 
lar IMF, increasing the fraction of high-mass stars by 



inhibiting the formation of low-mass stars. Theoreti- 
cal studies of the IMF all involve scaling from a Jeans 
mass, the mass at which a cloud's self-g ravity over- 
comes its support forces (e.g., iLarsonl [2005| ) . The sup- 
port can be either simple thermal pressure or large- 
scale turbulent motions . In th eir recent analytical study, 
iHennebelle fc ChabrieJ ()2008[ ) note that for typical ISM 
physical conditions and a reasonable prescription for the 
turbulence, the turbulent support is more important for 
the high-mass end of the spectrum, while the evolution 
of lower-mass condensations are governed by simple ther- 
mal support. 

The thermal Jeans mass can be written as Mj = 
l.lA/o(T/10i^)i Vr9°'^ (where Pig is the mass density in 
units of lO^^^gcm^'^) and yields 50 Mq for our median 
derived temperature (120 K) and density (10^'^ cm~'^). 
This may be indicating that the bulk of the material is 
indeed unlikely to participate meaningfully in any fur- 
ther SF in its present condition. Of course, the SF will 
occur in the densest and coolest regions, but they would 
likely be in approximate pressure equilibrium with the 
bulk of the gas. If we consider the lowest temperatures 
allowed by our likelihood analysis, T ~30 K, together 
with the highest pressures, log P ~ 6.7, then the density 
is n = 10^'^ cm""^, and the Jeans mass is Afj ~ 3Mo0 

A meaningful comparison with the Galaxy is hampered 
by the fact that our large-beam M82 observations are 
necessarily averaging over multiple SF regions, and will 
include gas in outflows as well as collapsing protostars 
themselves. Our approach is to examine the material 
around the Galaxy's most massive SF sites, since they are 
likely the best Galactic examples o f SF on large scales . 
We consider the sample studied bv lLeurini et al.l (|2007f ) 
in the millimeter and centimeter-band methanol transi- 
tions which are used to derive accurate temperatures and 
densities. We consider only the envelopes rather than the 
cores since the cores appear to be heated internally and 
are presumably already undergoing collapse, and in any 
c ase, the envelopes do minate the mass of these regions. 

iLeurini et al.l (|2007[ ) find temperatures ranging from 
11-36 K, and densities of 10^-10^ cm~^. Thus even these 
massive star formation sites are cooler on average than 
the dense gas in M82. The lower temperature is not 
surprising; again, in dense regions the gas temperature 
will approach the dus t temperature, whic h in the Galaxy 
ranges from 10-20 K (jParadis et al.il2009[ ). While the in- 
ferred thermal pressures in these star forming envelopes 
are comparable to those we find in M82, the lower tem- 
perature and higher density corresponds to a smaller typ- 
ical Jeans mass — values range from 0.3-2.8 M©, less than 
the minimum ~ 3M0 derived above for M82. 

More generally, the characteristic formed stellar mass 
scale M* is seen to scale as T^^^^, where Tmin is the min- 
imum temperature to which the gas can cool, and the 
exponent 7 ranges from 1.7 (o btained in numerical ex- 
periments, Jappscn et al. 2005) to 3.35 (via an analytic 
treatment. iLarson.l985i) . If we take the measured mini- 

^ We note for completeness that there is evidence for a (mass- 
independent) efficiency factor that relates the mass of a Jeans- 
unstable core to the mass of the actual star which forms from it, 
believed to be ~1.4-2 (HcnncbcUc & Cliabricr 2008), meaning that 
the resulting stellar masses are somewhat smaller than the Jeans 
mass estimates. 
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mum of ~30 K versus a conservative 20 K in the warm 
Galactic regions, this scahng suggests a factor of at least 
2-4 in for M82 relative to the Galaxy. Clearly, ac- 
curate estimates require theoretical study and more de- 
tailed knowledge of the local conditions at the star forma- 
tion sites, but if the gas we are tracing is indeed involved 
in star formation, then it likely produces a stellar IMF 
which is biased against low-mass stars relative to even 
the massive star formation sites in the Galaxy. 

Such a scenario is of course consistent with the reports 
of low-mass-dcficicnt stell ar pop ulations in M82 over the 
years (fRieke et all Il980l 119931: iForster Schreiber et all 
|2003[) . Moreover, an IMF biased against low- mass stars 
produces more luminosity per unit stellar mass than if 
the IMF is as observed in the Galaxy. Such a top- 
heavy or bottom-light IMF has been proposed to ex- 
plain an apparent discrepancy between the observed stel- 
lar mass buildup and the energy release history in the 
first h alf of the U niverse (z > 1) (jPerez-Gonzalez et al.l 
120081: iDavel I2008D . Given that the typical star form- 
ing galaxy in this epoch is now believed to be similar 
to the local LIRGs a nd ULIRGs (|Le Floc'h et al.l 120051 : 
iPapovich et al.ir2007l ). the conditions in M82 are likely 
more indicative of the historical average than those of 
the Galaxy. 

6. CONCLUSIONS 

We present a study of the dense molecular gas in the 
starburst nucleus of M82 based on 190-307 GHz spectra 
toward three positions obtained with the Z-Spec instru- 
ment. Z-Spec offers good sensitivity, accurate continuum 
measurement, and a uniform calibration for spectral lines 
across this band. We report fluxes for some 20 molecular 
transitions, many new detections. The measurements of 
the J = 3^2 transitions of HCO+, HON, and HNC, and 
the J = 4 — >■ 3 and J = 5 — 4 transitions of CS motivate 
an excitation and radiative transfer analysis in which all 
four species are simultaneously considered, incorporating 
all of their available published transitions. Our analysis 
constrains the physical conditions in the dense gas as well 
as the relative abundances among these species. We trace 
some 1.7-2.7x10® M© of gas with nu^ ~l-3x 10** cm'^, 
and find that it is warm: likely above 50 K and poten- 
tially as high as 500 K, a range which exceeds the level 
temperature of the transitions studied. The mass and 
temperature are thus comparable to that found for the 
warm component in the mid- J CO studies, but the higher 
density implies a thermal pressure of 1.5-4x 10^ K cm~^, 
about an order of magnitude higher then inferred from 
the mid- J CO transitions. 

In the framework of physical and chemical models, the 
line ratios among HCN, HCO+, and HNC indicate that 
the molecular gas is subject to both UV photons as well 
as a bulk heating mechanism other than X-rays. A sim- 
ilar conclusion is reached in considering the direct ob- 
served cooling in the CO lines up to J = 7 — >■ 6. Cosmic 
ray heating and dissipation of mechanical energy from 
the new star clusters are both potential heating sources 
for the molecular ISM in M82. This feedback has ren- 
dered much of the molecular ISM in the nucleus sterile to 
further SF. We briefly compare the dense molecular gas 
in M82 with star-forming sites in the Galaxy, concluding 
that if any of the material we are studying is involved 
in further SF, then the increased heating likely biases 



the stellar IMF against low-mass stars, relative to the 
Galaxy. Such a scenario may be more indicative of the 
typical SF environment in the Universe's history than 
the Galactic stellar IMF. 
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